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ABSTRACT 

Binary sequences have had application in communication 
systems for many years. Shift registers have been used in 
their generation, because of the ease and economy of their 
operation. For certain applications, nonlinear feedback 
functions are used by shift registers of span n to generate 
sequences of lengths up to 2 n . 

The sequences of maximum length 2 n and their generation 
are the subject of this thesis. In particular the ways 
of generating these sequences using nonlinear feedback 
shift registers and their correlation to linear feedback 
shift registers are described. Complexity is the term 
given to the length of the shortest linear feedback shift 
register generating a maximum length 2 n sequence. 

Games and Chan [Ref. 1] have given considerable study 
to the subject of complexity. Some of the problems they 
left are discussed further in this paper. It will be shown 
that the complexity of a de Bruijn sequence (S) is the same 
as the complexity of its reverse (r S) , complement (S’) , and 
its reverse complement (r S’) . Sequences (S) for which r S = 

S’ are termed RC sequences. It is shown that RC sequences 
exist for every odd n>_3 . In addition a lower bound will be 
established for the number of RC sequences occurring for 
each odd n. 
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INTRODUCTION 



I . 



A. APPLICATION 

de Bruijn sequences have been the subject of recurring 
interest since their discovery over a century ago. Recent 
developments in radar ranging [Ref. 2: Ch . 13], error correct- 
ing codes, secure or limited access code generators [Ref. 3: 
pp . 12-14], and mathematical modeling have made application 
of de Bruijn sequences because of their properties. The 
binary shift registers readily model digital computers with 
electronic states where "1" indicates the on state and ”0" 
indicates the off state. 



B. DEFINITION OF de BRUIJN SEQUENCES 

For a positive integer n, a de Bruijn sequence of span 

n is a complete binary cycle of length 2 n which is a sequence 

S = (s,, s 0 , •••, s } taken in circular order (s-, follows 
1 2 2 n 1 

s ) such that all possible n- tuples occur exactly once. 

2 n 

[Ref. 4: p. 120]. Complete binary cycles are often called 

de Bruijn sequences after the Dutch mathematician, N. de 

Bruijn, who in 1946 proved the existence of complete binary 

2 n_1 -n 

cycles having length 2 n and that they numbered 1 L 

de Bruijn was preceeded by Flye Saint-Marie a half century 

earlier in an obscure publication [Ref. 5]. As an example 

2 3 ‘ 1 -3 

for n = 3, Table 1.1 shows both (2 ) de Bruijn sequences 
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3 

of length eight (2 ) . Both of these sequences are examined 
to confirm that each of the eight 3-tuples occur exactly 
once. Notice that 2 elements of the sequence are repeated 
as signified by the bar ((TO) , which was necessary on each 
sequence to complete the last two 3-tuples. 



TABLE 1.1 

3-TUPLE COMPOSITION OF de BRUIJN SEQUENCES OF SPAN 3 



0001110100 

000 

001 

Oil 

111 

110 

101 

010 

100 



0001011100 

000 

001 

010 

101 

Oil 

111 

110 

100 



There are 2 n different possible starting points for 
each de Bruijn sequence; however, each cyclic permutation 
is considered equivalent. It is often a matter of conveni- 
ence to start each sequence in a canonical way with the all 
zero or all one n-tuple. This reduces confusion when working 
with a number of different de Bruijn sequences. 



C. PSEUDO- RANDOMNESS 

Since each binary n-tuple occurs exactly once on the 
cycle, the sequence models a uniform distribution with the 
probability of given n-tuple occurring to be exactly '*'/ 2 n * 
The sequence is not entirely random, however, since each 
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n-tuple (a^ a.^ ••• a fl ) has but two possible successors 

(0 a^ *** a n -i or 1 a i a 2 *** a n-l^’ ^ rua j n sequences 

satisfy the randomness properties of equidistribution of 
0's and l's and run lengths [Ref. 3: p. 10] which would be 
expected from the tossing of a fair coin. Thus, de Bruijn 
sequences have a pseudorandom property. 



D. K-ARY SEQUENCES 

In general it is possible to work with k possible states 

for each position of an n-tuple. Flye Saint-Marie in 1894, 

showed the existence and determined the number of complete 

,n-l . n-1 

cycles for n-tuples of k characters to be [(k-1)!] k 
[R ef. 5]. This report is restricted to the binary case 
(k = 2) due to its principal application in electronic com- 
munications and computers. 



E. DEFINITIONS 

For S = f s ^ S 2 ••• s^} where s^ e {0,1} for l£i<k the 
following are defined: 

1. Let W(S) denote the weight of S. 

k 

W(S) = S s- 
i=l 1 

2. Let P(S) denote the parity of S. 



P(S) = E s. (mod 2) 
i = l 1 
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3 . 



Let 1 (S) 

l(s) 



denote the length of S. 

= k where k is the number of positions in S. 
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II . SEQUENCE GENERATORS AND OPERATORS 



A. SHIFT REGISTERS 

In practice shift registers are used to generate de Bruijn 
sequences. An n-stage shift register has n memory registers 
(x^, x 7 , x n ) which shift their contents (0 or 1) to the 

next register upon command. For example, in Fig. 2.1, the 
contents of x^ will transfer to X 2 , x ^ to x^, etc., at the 
appropriate time with X n serving as the output. However, 
the contents of an n-stage shift register would empty in 
n shifts if no input was provided to x^ . 




Fig. 2.1 n-Stage Shift Register. 

B. FEEDBACK SHIFT REGISTERS 

Feedback networks are added to provide an input to the 

1st stage of the shift register. Note that in practice 

the output can be taken from any stage of the shift register. 

In fact, the contents of any stage is the same as the con- 

t h 

tents of any stage is the same as the contents of the n 1 " 



IS 



stage, merely shifted by a certain amount. In Fig. 2.2 the 
feedback function f is added to the n-stage shift register 
to generate non-trivial sequences. The coefficients 
a 0 , a^, a n are 0 or 1 with a Q = a n = 1. 




OUT 



Fig. 2.2 n-Stage Feedback Shift Register. 



1 . Linear Feedback Functions 

A feedback function f is classified as linear if 

the function f (x, , x 0 , •••, x ) is restricted to be of the 

1 2 n 

n 

form f = Z a 4 x- where the addition is modulo 2 addition. 
i = l 1 1 

The addition (mod 2) table is given in Table 2.1 below. 

Note that subtraction is equivalent to addition when opera- 
ting modulo 2. 
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TABLE 2.1 



ADDITION (MODULO 2) 



0 


0 


1 


0 


0 


1 


1 


1 


0 



2 . Nonlinear Feedback Functions 

If the feedback function f(x.,, x ) utilizes 

1 n 

multiplication (mod 2) and addition (mod 2) then the feed- 
back function is nonlinear. (We also say that the linear 
feedback functions are vacuously nonlinear feedback functions-) 
Modulo 2 multiplication is identical to multiplication in 
the integers with the restriction of the alphabet used 
being the set (0, 1} as described in Table 2.2. 

TABLE 2.2 

MULTIPLICATION (MODULO 2) 



• 


0 


1 


0 


0 


0 


1 


0 


1 



3 . Comparison of Linear and Nonlinear Feedback Functions 
Considerable research has been devoted to the study 
of linear feedback functions because of their ease of analy- 
sis. Nonlinear feedback functions are much more complicated. 
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However, the class of sequences that can be generated by 

nonlinear functions is much greater than for linear functions. 

The number of binary maximum length (2 n -l) linear 

sequences that can be generated by a n-stage linear feed- 

2 n 

back shift register (LFSR) is approximately — [Ref. 7]. 

The reason that the maximum length linear sequences only 

have length 2 n -l is that the zero state is a fixed point 

under the linear function. This poses no problem since the 

maximal length linear sequences are readily made into 

de Bruijn sequences having length 2 n simply by adding a 

zero to the (n-1) -tuple of zeros. The feedback function 

involved then becomes a nonlinear function. 

By comparison the number of binary nonlinear full 

length sequences (2 n ) that can be generated by a nonlinear 

2 n-l 

feedback shift register (NFSR) is exactly 1 L . Below, 

Table 2.3 compares these numbers for 3<n<8 . The number of 
maximum length sequences will be of interest later in 
Section III, where the concept of complexity in generating 
these sequences will be developed. Various algorithms are 
listed in a survey by Fredricksen [Ref. 5] for generating 
de Bruijn sequences. 

C. SUCCESSOR AND PREDECESSOR STATES 

The contents of an n-stage shift register at a specific 
time is usually referred to as its "state". As mentioned 
in Section I when discussing the pseudorandomness of de Bruijn 
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TABLE 2.3 



NUMBER OF MAXIMUM LENGTH SEQUENCES 



n 


Linear 


Nonlinear 


3 


2 


2 


4 


4 


16 


5 


6 


2048 


6 


6 


67,108,864 


7 


18 


17 

1.4 x 10 17 


8 


30 


2.7 x 10 j6 



sequences each state has two possible successor (conjugate) 
states depending on whether the feedback function generates 
a "0" or "1" as the input to register x^. Likewise each 
state has two predecessor (companion) states depending on 
whether the previous output was a "0" or "1". The adjacency 
quadruple in Fig. 2.3 centers around then (n-l)-tuple x 9 **»x n . 
That is, each input state has two possible successors depend- 
ing on whether the feedback is 0 or 1, and each output state 

has two possible predecessors depending on whether the pre- 

th 

vious output from the n stage was 0 or 1 . 

In Table 2.4 we show an example for n = 3 of the possible 
successors and predecessors for each state. 
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TABLE 2.4 

CONJUGATE AND COMPANION STATES (n=3) 



Predecessors 


State 


Successor 


000-001 


000 


000-100 


010-011 


001 


000-100 


100-101 


010 


001-101 


110-111 


Oil 


001-101 


000-001 


100 


010-110 


010-011 


101 


010-110 


100-101 


110 


011-111 


110-111 

t T 


111 


011-111 
t T 


Output 




Input 



D. de BRUIJN DIAGRAMS 

The de Bruijn diagram compactly contains all the above 
information. The diagram contains 2 n vertices corresponding 
to the 2 n states and two directed edges from each state to 
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the possible pair of successor vertices. As a result two 
arrows exit every vertex leading to conjugate states, and 
two arrows enter every vertex arriving from companion states. 
For an example consider the de Bruijn diagram G n for n = 3 
in Fig. 2.4. An Eulerian path through G n is defined as a 




Fig. 2.4 Graph G^. 
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path which visits each edge exactly once. A path which 
visits every vertex of G n is called a Hamiltonian path. A 
Hamiltonian path in G n is a de Bruijn sequence of span n. 

If the edges of G n are labeled with an n + 1-tuple defined 
by the labels on the predecessor and successor states, then 
the Eulerian paths in G n correspond to Hamiltonian paths-- 
and de Bruijn sequences in the graph G n+ ^. The Hamiltonian 
path in Fig. 2.5 is a subgraph of G~ in Fig. 2.4 with exactly 
one edge emanating from each vertex. 

E. OPERATORS 

For S = s, s~ ••• s . s where s. e {0,1} for l<i<n the 
12 n-ln l — — 

following operators are defined. 

1 . Identity Operator 

The identity operator (e) operating on S denoted 

eS = s, s~ • • • s . s =S. 

1 2 n-1 n 

2 . Reverse Operator 

The reverse operator (r) operating on S is denoted 

r S, where r S = s s , ••• s 9 s,. The reverse operator 

n n- 1 21 r 

preserves the weight, parity and length of S. Note also 

2 

that (r) = e. 

3 . Complement Operator 

The complement (dual) operator on S is denoted d S 

or S, where S = s, s 0 ••• s’ , s’ such that s’. = s. @1 

12 n-ln ii 

(mod 2) for l<i<n. 
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Fig. 2.5 A Hamiltonian Path through G,. 
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s n then 



Lemma 1: If W (S) = k for S = s ^ s ? ••• 

W (S) = n-k. 

Proof: W (S) = k implies k nonzero elements z S 

and n-k zero elements z S. 
Complementing n-k zero elements produces 
n-k nonzero elements z S. 



I . W (S) = n-k 

Q . E . D . 

Lemma 2: If P (S) = a where a z {0,1} and S = {s, s~ ••• s }, 

1 2 n 

then P (S’) is a iff n is even, 
n 

Proof: P (S) = Z s. (mod 2) = a 

i = l 1 



n _ n 

P (S’) = Z s • (mod 2) = Z s - @ 1 = a @ n 
i=l 1 i=l 1 



therefore P (S’) = P (n-a) 



Q.E.D. 

2 

Also note (d) = e, and that the reverse and comple- 

ment operators are commutative, i.e. rd = dr. 

4 . Reverse Complement Operator 

The reverse complement operator on S is denoted 
r S^ = s” n ? n _ 1 ••• s, . Weight and parity are effected 
by the reverse complement operator exactly as with the 
complementation operator. 

These operators apply to sequences as well, since 
sequences are composed of n-tuples. Table 2.6 shows the 
effect of operators on de Bruijn sequences of span 3 and 4 
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where r S, S, and r S are all de Bruijn sequences, as they 
always are if S is a de Bruijn sequence. 



TABLE 2.5 

de BRUIJN SEQUENCE OPERATORS 



n = 3 

S = 00011101 = 

rS = 10111000 = 

S = 11100010 = 

rS = 01000111 = 



(Canonical Form) 
00011101 
00010111 
00010111 
00011101 



n = 4 

S = 0000111101001011 
rS = 1101001011110000 
S = 1111000010110100 
rS = 0010110100001111 



(Canonical Form) 
= 0000111101001011 
= 0000110100101111 
= 0000101101001111 
= 0000111100101101 



If as in the case for n = 3 in Table 2.5 a sequence 
S is equivalent to its reverse complement, then S is termed 
a RC sequence. This applies to both de Bruijn and non- 
de Bruijn sequences. The following theorem appears in a 
paper by Etzion and Lempel [Ref. 6]. 

Theorem 1: A sequence S is a RC sequence if l (S) is 

even and S = {X, rX} for some X where X is a binary string. 
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F. LINEAR RECURSION FORMULA 



de Bruijn sequences generated by a nonlinear function 
can also be generated by a linear function f (x^ ••• x n ) = x Q 
having the linear recursion formula 
n 

Z a . x • = 0 . 

j =0 J J 



Here each a^ determines a tap to the j 



, th 



stage of the 



register. The linear generator polynomial must be 
(x + 1) for some positive integer k, due to the periods of 
de Bruijn sequences (2 n ) and results from theorems given 
by Golomb [Ref. 7: pp. 27-43]. 

1 . Pascal Triangle (Mod 2) 

The values for each a^ are summarized in the Pascal 
triangle (mod 2) Table 2.6 for various values of n. They 
are the coefficients generated by the binomial expansion of 
(x + 1)^ mod 2. 



TABLE 2.6 

PASCAL TRIANGLE (MOD 2) 

k = 0 1 

1 11 

2 10 1 

3 1111 

4 1 0 0 0 1 

5 110011 (x+1) 5 

a a. a, a^ a. a c 

o 1 2 o 4 5 
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2 . Pascal Triangle Properties 

(a) The a^ sequence is palindromic, that is 

= a n = !, a i = a n _ ]_ > a 2 = a n-2 



etc . 



(b) The number of nonzero coefficients is even, 
n 

that is E a. = 0 (mod 2). 

j=0 J 
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Ill . COMPLEXITY OF de BRUIJN SEQUENCES 



A. DEFINITION OF COMPLEXITY 

The term complexity indicates a general measure of the 
predictability of a sequence. Various definitions of com- 
plexity are in use, but this paper will use the one given 
by Chan and Games [Ref. 1]. Their definition of complexity 
is the length of the shortest LFSR required to generate a 
sequence. They apply their results primarily to de Bruijn 
sequences. This definition is also the one used by Herlestom 
[Ref. 8] to look at two different shift register cycle gen- 
erators. For a de Bruijn sequence S, let C(S) denote the 
complexity of S. 

Previous results by Chan and Games [Ref. 1] establish 
lower and upper bounds of 2 n ^+n and 2 n ^ respectively for 
the complexity of a de Bruijn sequence. Thus, the length 
of a LFSR is almost as long the sequence as opposed to a 
n-stage NFSR which generates the sequence. Table 3.1 serves 
as a ready reference for the upper and lower limits of 
complexity for 3<n£ 8 . The upper bound is known to be attained 
for all n and the lower bound for all n<6 . 

B. COMPLEXITY ALGORITHM 

A fast method for determining the complexity of a de 
Bruijn sequence was developed by Games and Chan [Ref. 9] 
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TABLE 3.1 



Span 

(n) 

3 

4 

5 

6 

7 

8 



LIMITS OF COMPLEXITY 



Lower Bound 
j2 n ~ 1 *n) 

7 

12 

21 

38 

71 

136 



Upper Bound 
= (2 n -l) 

7 

15 

31 

63 

127 

255 



and is represented by the flow chart in Fig. 3.1. Let S be 



a de Bruijn sequence then &(S) = 2 n . Further let S = A:B = D 



where A = {a^ a 2 



a , } and B = {b, b 0 
2 n ' 1 1 2 



b , } then 
2 n ~ 1 



D = {a x a 2 



a . b i b, 
2 n ” 1 1 2 



b , } . To serve as a guide 
2 n- 1 ° 



an example for n = 5 is now presented. Let the de Bruijn 
sequence S = {11111000001000110101100101001110} having length 



i = 32. 

A 1111100000100011 a = 16 

B 0101100101001110 C = 0 

A+B 1010000101101101 C = 16 

A 10100001 l = 8 

B 01101101 C = 16 

A+B 11001100 C = 24 
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A 1100 
B 1100 
A+B 0000 



l = 4 Since A€>B = 0 

C = 24 

C = 24 C is unchanged 



A 11 
B 00 
A+B 11 



l = 2 
C = 24 
C = 26 



A 1 
B 1 
A+B 0 



l = 1 Add 1 since l = 1 

C = 26 fj A = 1 

C = 26 + 1 = 27 



C. COMPLEXITY DISTRIBUTION 



The nature of the distribution of complexity of de Bruijn 
sequences is the primary interest of this paper. Let 
a (c,n) denote the number of de Bruijn sequences having com- 
plexity c and span n. The complexity distribution of de 
Bruijn sequences of span n for 3<n<_6 are listed in Table 

3.2. Unfortunately the large number of de Bruijn sequences 
5 7 

for n>_7 (2 ) does not allow an exhaustive examination of 

their complexities. 

Examination of this data led Chan and Games to conjecture 
that for n>_3 a (c,n) is congruent to 0 modulo 4. 



D. EFFECTS OF OPERATORS ON COMPLEXITY 

Chan et al. [Ref. 1] proved that a(c,n) for n>_3 was con- 
gruent to 0 modulo 2. By showing that for each de Bruijn 



30 







* 



C (S) =c 



(^ STOP^ ^ 



Fig. 3.1 Flowchart for Complexity Algorithm. 
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TABLE 3.2 

COMPLEXITY DISTRIBUTION 



c 


n = 4 
a(c,4) 


c 


a(c,S) 


n = 5 

c 


12 


4 


21 


*• 


27 


13 


0 


22 


0 


28 


14 


4 


23 


12 


29 


15 


8 


24 


20 


30 






25 


32 


31 






26 


36 





n = 6 



c 


a(c,6) 


c 


a (c , 6) 


38 


448 


51 


8704 


39 


o • 


52 


18096 


40 


32 


53 


34224 


41 


96 


54 


67700 


42 


160 


55 


126592 


43 


80 


56 


259320 


44 


432 


57 


519752 


45 


288 


58 


1041252 


46 


896 


59 


2090716 


47 


1168 


60 


4162352 


48 


2772 


61 


8342176 


49 


2352 


62 


16692832 


50 


5224 


63 


33731200 



(c , 5) 

64 

180 

224 

448 

1024 
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sequence having complexity c, there exists another different 
de Bruijn sequence S also having complexity c. Consequently 
they showed that the complement operator does not change the 
complexity of a de Bruijn sequence. Using a different 
approach, the next theorem will show that the complexity of 
a de Bruijn sequence is also the complexity of the de Bruijn 
sequences defined by the reverse operator and by the reverse 
complement operator as well. 

Lemma 3 . For S = {s-^ S 2 ••• s n _^} , a de Bruijn sequence, 
c 2 

C(S) = c iff Z a- s. = 0 for everv s. where each a. comes 

th i=0 1 1 1 

from the c row of the Pascal triangle (mod 2) and c is the 

smallest integer such that the above recursion holds. 

Proof . S is de Bruijn with C(S) = c iff s^ = f(s^ + ^*»* 

s- „) for every s. e S, which is equivalent to s • = 
i+c' ' 1 n 1 

c 

Z a- s i+ . from results by Golomb [Ref. 5: pp . 27-43] 
j =1 J 1 J c 

mentioned earlier. In turn s. = Z a- s- • = 0 from proper- 

1 j =0 ^ ^ 

ties of the Pascal triangle given in Section II. Q.E.D. 

This lemma is used in the following theorem to show that 
complexity is preserved by various operators on de Bruijn 
sequences . 

Theorem 2. If S is a de Bruijn sequence having complexity 
C, then C(S) = 0(3") = C(rS) = C(r3) = c. 

Proof . Part I (reverse operator) 

Let s^ *** s i+ c be an arbitrary (c+1) long segment 

of S which satisfies 
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0 



(1) Z a. s 



j-o 



J i + J = 



by Lemma 3. Since the a^ sequence is palindromic as shown 

earlier, s^ +c ••* s^ + ^ a (c+1) long segment of rS satisfies 

equation (1). Thus C (rS) = c. 

Part II (complement operator) substituting ?^ ••• sL +c 

the c+1 long segment of S’, into equation (1) , the results are 

a ?-+ a, ?-,+•••+ a , ?. , + a ?. 

oi 1 l+l c-1 i+c-1 c l+C = 



a o ( - s i + ^ + a l ^ s i + l + 1 ) + **• 



+ a 



c-1 



(s i-c-l + 1) 



+ a 



c c 

(s. „ +1) = Z a. s. . + £ a. (mod 2). Since there are 

1+C j =0 3 1 + 3 j=0 3 

c 

an even number of nonzero a.. s, E a. = 0 and equation (1) 

J j=0 3 

is satisfied. Therefore, by Lemma 3 C(ST) = c. 

Part III (reverse complement operator) , by Parts I and 
II, it follows immediately that C (rS”) = c. Thus C (S) = 

C (rS) - C (?) = C (rS) 



Q.E.D. 
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IV. TRUTH TABLES 



Linear and nonlinear feedback functions were used in 

conjunction with n-stage shift registers to generate de Bruijn 

sequences in Section II. The input was generated as some 

function f of the current state x, x 0 x„ , that is 

1 z n ’ 

(1) X, = f (Ij X; X n ) 

or equivalently writing f as in (2) when specializing for 
cycles only functions. 

(2) x o = x n * 8 Cx l ••• Vd- 

Extensive use is made of equation (2) since results for 

cycles only are desired. 

A. FULL TRUTH TABLE 

The full truth table is a useful way of listing the 
functional value Xq from equation (1) for each possible 
state of the shift register. The functional values in Table 
4.1 is an example of a full truth table for a de Bruijn 
sequence (00011101) . 

Examining the full truth table it becomes apparent that 
the information in the lower half (below the dashed line) 
is redundant, since the values of f (x^ X 2 x^) for x^ = 1 
in the lower half are merely the complements of the values 
of f (x.^ x 2 x 3 ) for x 3 = 0 in the upper half for the identical 
values of x^ x 2> This suggests that the upper half (or the 
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TABLE 4.1 



FULL TRUTH TABLE OF (00011101) 



Current 
(x 3 x 2 


State 

x x ) 


Functional Value 
x 0 = f (x 1 x 2 X,) 


0 


0 


0 


1 


0 


0 


1 


1 


0 


1 


0 


0 


0 


1 


1 


1 


1 


0 


0 


0 


1 


0 


1 


0 


1 


1 


0 


1 


1 


1 


1 


0 



lower half) of the truth table can display all the informa- 
tion for a function generating only cycles. 

B. HALF TRUTH TABLE 

Focusing on equation (2) , since x n is always 1 in the 
lower half this can be viewed as complementing the functional 
values for g (x^ ••• x in the upper half. In this manner 

the half truth table in Table 4.2 is constructed for the same 
de Bruijn sequence (00011101) to aid comparison with Table 
4.1. From this point on the term truth table will always 
mean the upper half truth table with x n = 0 , unless noted 
otherwise . 
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TABLE 4.2 



HALF TRUTH TABLE OF (00011101) 



X 2 


X 1 


g (x x x 2 ) 


0 


0 


1 


0 


1 


1 


1 


0 


0 


1 


1 


1 



The string of functional values in the truth table is 
defined to be the generator G = (g^ g^ ••• g n _ x } where 
the subscripts of g are the decimal equivalent of 



k n-l 



Both the weight and parity of generator G 



will be of interest as the investigation of this paper con- 
tinues. Whenever the parity or weight of the truth table 
is mentioned in this paper, it is referring to the generator 
of that truth table. 



C. ZERO TRUTH TABLE 

The zero truth table is defined as the truth table in 
which all functional values are zero. The zero truth table 
for n = 3 contained in Table 4.3 is basically the output of 
a circulating register which generates pure cycles. 

Pure cycles are defined as cyclic permutations of the 
original state. The zero truth table in Table 4.3 generates 
the four pure cycles shown in Fig. 4.1. The number of pure 
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TABLE 4.3 



ZERO TRUTH TABLE FOR n = 3 

x 2 x 1 g (x 1 x 2 ) 

0 0 0 

0 1 0 

10 0 

11 0 

cycles Z has been shown [Ref. 7: p. 120] to be 

Z n = 1 Z 4> (d) 2 n/d . 
n d| n 

Eulers totient function <j) (d) is the number of fractions 

of the form a /d where l£a<d and a /d is in lowest terms. The 
summation is over all positive ‘integers d which divide n, 
denoted (d|n). It is further shown [Ref. 5] that Z n -1 is 
the minimum possible weight truth table that can generate a 
de Bruijn sequence. 

In the next section, the proof of theorem 6 gives an 
example of how pure cycles can be joined in a way to create 
a de Bruijn sequence that has a minimum weight truth table. 
Table 4.5 gives a listing of Z n and Z n ~l for l£n<_7 . 

D. ONES TRUTH TABLE 

The ones truth table assigns a functional value of 1 to 
every entry of g (x^ ••• x n .j) • The ones truth table depicts 
the output of a complementing shift register. The input is 
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Fig. 4.1 Pure Cycles (n=3) . 
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always the complement of the output, i.e. f (x^ ••• x n ) = 

1 + x . The ones truth table for n = 3 is shown in Table 
n 

4.4. The complementing cycles from Table 4.4 are depicted 
in Fig . 4.2. 



TABLE 4.4 

ONES TRUTH TABLE (n = 3) 




The number of complementing cycles Z n .[Ref. 7] is 



n 



= Z - 1 E 4> (2d) 2 
— y 2n 2d|n 



7X 



Since summation for Z is over all 2d (even) numbers which 

n v J 



divide n, for n odd Z 



n 



= Z. 



_n 

2 



Accordingly Fredricksen [Ref. 5] shows truth tables of max 
weight that can generate a de Bruijn sequence is 



,n- 1 



n 



+ 1 . 



The values of Z * and 2 n ~^ - Z +1 are given in Table 4.5 

n n 

for 1 <tk 7 . Also included are the number of de Bruijn se- 
quences having truth tables of maximum or minimum weight 
from a listing by Fredricksen [Ref. 5]. 
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TABLE 4.5 



PURE CYCLE COMPILATION 



n 



Zn 



V 1 



# de Bruijn 
Sequences 
with Min 
Wt Truth Table 



V 



2 n_1 -Z * + l 
n 



# de Bruijn 
Sequences 
with Max 
Wt Truth Table 



1 

2 

3 



2 

3 

4 



1 

2 

3 



1 

1 

2 



1 

1 

2 



1 

2 

3 



1 

1 

2 



4 6 5 



12 



2 7 



3 



5 8 7 

6 14 13 

7 20 19 



2 6 *3 2 
2 14 • 3 4 • 5 2 
2 28 *3 5 -5 3 *13 



4 

6 

10 



13 

27 

55 



2 



6 



2 



14 




• 3 
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Fig. 4.2 Complementing Cycles (n=3) 

E. EFFECTS OF CHANGING THE TRUTH TABLE 

If the ones truth table in Table 4.4 is changed 
g (0,1) = 0 vice 1, the truth table of Table 4.6 is 



so that 
produced . 
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TABLE 4.6 



TRUTH TABLE (TABLE 4.4 MODIFIED) 



X 2 


X 1 


g (x x x 2 ) 


0 


0 


1 


0 


1 


0 


1 


0 


1 


1 


1 


1 



This change in the truth table effects the cycles of Fig. 

4.2 by joining them as depicted in Fig. 4.3. Indeed the 
change produces a Hamiltonian path for the de Bruijn sequence 
( 00010111 ) . 

Suppose the truth table of Table 4.6 is now changed so 
that g (00) = 0 vice 1. This produces the truth table in 
Table 4.7, and its cyclic composition in Fig. 4.4. 

TABLE 4.7 

TRUTH TABLE (TABLE 4.6 MODIFIED) 



x 2 


X 1 


g (x x x 2 ) 


0 


0 


0 


0 


1 


0 


1 


0 


1 


1 


1 


1 
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g. 4.3 Cyclic Composition of Table 4.4. 
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Fig. 4.4 Cyclic Composition of Table 4.7. 
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This change in the truth table causes the cycle to split 
into two cycles. Suppose (x ••• x^) is on element on upper 
of two cycles as shown in Fig. 4.5 and that its successor is 




Fig. 4.5 Two Cycle Graph. 



X , • • • x, 1 as f (x. ••• x ) = 1. If the truth table is 
n-1 1 1 n ' 

now changed so that f (x. ••• x ) =0, then x , ••• x, 0 
° In n-l i 

is the new successor for x n ••• x^ . Retaining the cycles 

only condition f fx, ••• x , x )must now be 1 . So the 
' 1 n-l rc 

successor of x x , • • • x, is now x , • • • x, 1 . One of 

n n-l 1 n-l l 

two things can happen how depending on whether x n ••• x-^ 0 

is on the upper or lower cycle. 
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First if ( x n _i **• x^O) is on the upper cycle, then that 
cycle is split as shown in Fig. 4.6 forming three cycles. 





Fig. 4.6 Cycle Splitting. 

Otherwise if (x ^ * * * x^O) is on the lower cycle, then 
those two cycles are joined as shown in Fig. 4.7. 

Clearly the parity of the cycles changes for each change 
in the truth table. The following theorem from Golomb 
[Ref. 7: p. 122] shows that the parity of the number of cycles 
and the parity of the truth table generator are equal. 

Theorem 3 . For n>2, the parity of the truth table is 
equal to the parity of the number of cycles of the truth 
table . 
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Since a de Bruijn sequence is one cycle by definition, 
the following corollary of theorem 3 is given. 

Corollary . If a sequence is de Bruijn, then the parity 
of the TT of that sequence is odd. 







Fig. 4.7 Cycle Joining. 

F. SPECIAL NUMBERS 

Analyzing the effects of de Bruijn sequence operators 
on their truth tables, certain binary numbers are found to 
have special properties which can be exploited in later 
analysis . 
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1 . Type RC Numbers 
If X = x, •« 



n 



X 1 = X 1 



• x n is a binary number, such that 

x„ then X is called a RC number, 
n n 



Lemma 4 . For n even there are exactly 2“ type RC numbers. 

Proof . For n even, if X is an RC number, then 

:, • • • x = x ••• x, bv definition. Therefore, 

1 n n 1 

= x , x- = x , , etc. Thus x,***x can be written as 
1 n Z n-1 1 n 



x, x~* • *x x • 
12 n n 

2 I 

: t 
n 

,7 



n 



X 2 ^-^ . Since j positions can be filled in 



either of^two ways, the total number of RC numbers for n 



even is 2‘ 

Q.E.D. 

There are no type RC numbers for n odd, since this 
would require t>hat x n+ ^ = x n+ ^ which is impossible. 

2 . Type R Number^ 

If X = x i*** x n is a binary number such that 

x, *»*x = x •••x 1 , then it is called a type R number. 

1 n n 1 IL 

2 

Lemma 5 . For n even there are + |xactly 2 type R numbers, 

— 9 

and for n odd there are exactly 2 type R numbers. 

Proof . The proof for n even parallels that of Lemma 4 
and will not be repeated. For n odd, since x n+ ^ = x n+ ^ 
there are n+1 positions to be filled in either"g|-^two^ 



ways. Therefore, for n odd there are exactly 2 
numbers . 



type R 



Q.E.D. 



49 



A useful way of organizing special numbers and other 
truth table values to assist in evaluating effects on the 
generators for various operators is shown in Table 4.8 for 
n = 5 and n = 6. The numbers are the decimal equivalent of 
the binary (n-1) -tuple, from the feedback function 

x 0 * x n + 8 (X 1 "• x n-l J 
which also corresponds to the truth table state. 

For example for n = 5, if f (2) = 1 for some sequence 
S then: f (11) = 1 for r S, f (4) * 1 for r S, and f (13) = 1 
for S. In general, for some operator 8 and some number a, 
if f (a) = 1 for some sequence S, then f (8a) = 1 for 8 S. 
Special numbers have properties in addition to this. If a 
is a type RC number and f (a) = 1 for S, then f (a)= f (r a) = 
1 for r S. Similarly if y is a type R number and f (y) = 1 
for S, then f (y) = f (ry) = 1 for r S. Thus, if we want to 
shown S is a sequence such that S = r S’, then if f (a) = 1 
in S, then f (r a) = 1 in S also. 

G. GENERATOR ANALYSIS 

Analysis of the generators for the de Bruijn sequences 
S, r S, S, 5 r S in Table 4.9 appears to indicate the gener- 
ator G of sequence S denoted G (S) is the reverse of the 
generator G (S) . Notice also that G (r S) = r [G (r S ) ] 
for S. Although Table 4.9 is an example of a single sequence 
S, it is not an isolated incident; in general the result 
G (S) = r [G (S) ] holds and is proven in the following 
theorem. 
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TABLE 4.8 



TRUTH TABLE ANALYSIS CHART 









rS 




s 


S 


rS 




s 


S 


rS 


rS 




rS 


rS 


S 


S 




rS 


rS 


S' 


S 




0 


15 


0 


15 


Truth table states 


1 


14 


8 


7 


Type R numbers 


2 


13 


4 


11 


(0,6,9,15) 


3 


12 


12 


3 


Type RC numbers 


5 


10 


10 


5 










(3,5,10,12) 


6 


9 


6 


9 




0 


31 


0 


31 


**- Truth table states 


1 


30 


16 


15 


Type R numbers 


2 


29 


8 


23 


(0,4,10,31) 


3 


28 


24 


7 


Type RC numbers 


4 


27 


4 


27 


(None) 


5 


26 


20 


11 




6 


25 


12 


19 




9 


22 


18 


13 




10 


21 


10 


21 




14 


17 


14 


17 





(n=5) 



(n=6) 
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TABLE 4.9 



GENERATOR ANALYSIS FOR de BRUIJN SEQUENCES 

(r S f S) 



S 

S 

rS 

rS 



11111001000111011010011000001010 

11111010100000110111000100101100 

11111010100000110010110111000100 

11111001101001000111011000001010 



X i 

n-1 


• • • 


,x 1 G (S) 


G (?) 


G (rS) 


G (r?) 




0 


1 


1 


1 


1 




1 


0 


1 


1 


0 




2 


1 


0 


0 


1 


□ 


E 


1 


0 


0 


1 


Type RC 


4 


0 


1 


1 


0 




E 


0 


1 


1 


0 


o 


© 


0 


1 


0 


1 


Type R 


7 


0 


1 


1 


0 


/other than\ 


8 


1 


0 


0 


1 


Vo § is y 


© 


1 


0 


1 


0 




LilJ 


1 


0 


0 


1 




ii 


1 


0 


0 


1 


1 


1 12 


0 


1 


1 


0 




13 


0 


1 


1 


0 




14 


1 


0 


0 


1 




15 


1 


1 


1 


1 
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Theorem 4 : If S is a de Bruijn sequence with generator 



G (S) , then S is a de Bruijn sequence with generator 
r [G CS) ] . 

Proof : Let G (S) = {g^, g 2 » g n _^} then the full 

truth table for S is given in Table 4.10. 



TABLE 4.10 

FULL TRUTH TABLE FOR S 



x n 


x n- 1 


X_ 9 * * • 

n - 1 


X 1 


f (x^ , •••, x n ) 


0 


0 


0 • • • 


0 


§1 


0 


0 


0 • • • 


1 


S2 


• 

• 

• 

0 


1 


1 • • • 


• 

• 

• 

1 


• 

• 

• 

V' 1 


1 


0 


0 • • • 


0 


*1 


1 


0 


0 • • • 


1 


g 2 


• 

• 

• 

1 


1 


1 • • • 


1 


• 

• 

• 

V" 1 ■ 



If x n , x nl , x lf f (x 1 ••• x n ) is a (n+1) long string 
in S, then x n , x n _ 1 , •••, Xj , f ^ ~ 2 * * • x n ) is a (n+1) 
long string in S’. 

Therefore we can transform each line of the full truth 
table for S by taking its complement to get Table 4.11. 
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TABLE 4.11 



FULL TRUTH TABLE FOR S (NON STANDARD] 




Rearranging the rows from top to bottom, the TT for S is 
placed in standard form as depicted in the Table 4.12. 

Further examination of the generators in Table 4.9 shows 
the structure of G (S) to be remarkably similar to the 
structure of G (r S) , with the exception that the functional 
values of type R numbers (6 5 9) are transposed. If, however, 
these values should be the same then clearly G (S’) = G (r S) 
or S = r S’. Indeed, this occurs twice for n = 3 and 64 times 
for n = 5, and will be shown to happen for all odd n>_3 . 

This cannot happen for n even, because there are no type 
RC numbers in the truth table. 
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TABLE 4.12 



Clearly 




(S) is r [G (S) ] 

Q . E . D . 
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V. RESULTS 



Examining Table 3.2 it can be quickly determined that the 
number of de Bruijn sequences of given complexity for 3_<n<6 
is 0 (mod 4). Games and Chan [Ref. 1] made this conjecture, 
but were able to show only that the distribution a (c, n) 
was 0 mod 2. Their result followed by showing that for a 
given de Bruijn sequence S of complexity c, there existed 
a unique de Bruijn sequence S also having complexity c. 

In Section III it was shown that reverse and reverse dual 
operators preserve the complexity of de Bruijn sequences as 
well . 

If each of these new operators produced unique de Bruijn 
sequences, then it would be easy to show that the numbers 
a (c, n) = 0 (mod 4) since for every de Bruijn sequence S 
there would also exist S’, r S and r S’ all de Bruijn and 
unique having complexity c. As mentioned in Section IV, the 
sequences r ST and r S are not always distinct from the 
sequences S § S’ respectively; however, for n even the next 
theorem shows that all four of these sequences are distinct. 
In addition it will be shown that RC sequences occur for 
all n odd, and a lower bound for the number of RC sequences 
will be established for all n and investigated for n = 7. 
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A. a (c, n) E 0 MODULO 4 FOR n EVEN 

Etzion and Lempel [Ref. 6] give a proof that a (c , n) 

is congruent to 0 mod 4 for c n even. In the proof for n 

even given here a different approach is taken. 

Theorem 5 . For even n>_3 , the numbers a (c, n) e 0 (mod 4) 

for a de Bruijn sequence S of span n and complexity c. 

Proof . From theorem 1, it is known that for a de Bruijn 

sequence S of complexity c, that S’, r S § r S’ also have 

complexity c. Games and Chan showed the uniqueness of S’, 

so all that remains is to show that r S is distinct from 

S, or equivalently that r S is distinct from S for n even. 

Suppose on the contrary that S = r S". Then consider 

an n-tuple Oa , ••• a, in the truth table such that 
r n- 1 1 

f (a^ ••• 0) = 1. Then 0a n ••• a^ 1 is a n+1 long 

string in S, and o a^ ••• a ^ 1 is a n+1 long string in r S’. 
If S = r S 

then 0 a^ ••• a^ _^1 is a n+1 long string in S, 
which implies 

that 1 a^ ••• a n _i 0 is a n+1 long string in S, and 

therefore f (a^ ••• a n _^ 0) = 1. This must hold for every 

x-^ ••• x n _i 0 such that f (x-^ ••• x n 0) = 1. Pairing 

these vectors, a, • • • a , and a . • • • a. , two at a time 

'1 n-1 n- 1 1 

results in an even number of l's in the truth table. So 

there must be a vector which pairs with itself as the parity 

of ones is odd. But then x . • • • x. = x. • • • x . is a 

n-1 1 1 n-1 
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contradiction since x f x J _, for n even. Therefore for 

n+1 n+1 

2 2 

n even S ^ r S’ and r S f S , so r S and r are distinct 
and a (c, n) e 0 mod 4. 

Q . E . D . 

B. RC SEQUENCES EXIST FOR ALL ODD n>3 

The results a(c,n) = 0 (mod 4) for all n are not obtained 
as hoped for. For n = 3 and n = 5 RC sequences are kno\vn 
to exist, but it is not known if this is true for all odd 
n>_3 or just some. The following theorem shows that RC 
sequences exist for all odd n>_3. Recall that RC sequences 
are those sequences S for which S = r S or equivalent by 
S = r S since r (S) = r (r S’) yields S = r S. 

Lemma 6 . The sequence formed by joining all pure 
cycles of n-tuples is de Bruijn. 

Proof . The least elements of each pure cycle for a 
given n are easily composed and arranged in order by weight 
according to the algorithm in Appendix A. Since the pure 
cycles contain all n-tuples exactly once, all that remains 
is to show how each pure cycle can be interconnected, thus 
forming a de Bruijn sequence. For each cycle of wt a>_l 
there exists a least element "2X+1" (necessarily odd, if 
it were even (2X) then "X" is on the same cycle and smaller) 
whose predecessor M X" has wt (a-1) . Since the pure cycles 
are ordered by weight, it is clear that the least element 
of each pure cycle of wt a^l can be joined to a cycle of 
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wt (a-1) by changing the truth table for its predecessor, 
i.e. g(X) becomes 1. In this way all the pure cycles are 
now interconnected, and the sequence is de Bruijn. 

Q.E.D. 

Theorem 6 . For all odd n>_3, there exist RC de Bruijn 
sequences . 

Proof . An RC de Bruijn sequence can always be composed 

in the following way for n>_3. First construct a light cycle 

n 1 

(LC) for those pure cycles of wt a<2-j— identical to the 

method used in Lemma 6. Note that when the cycles of 

wt 2^- are joined to cycles of smaller weight the truth 

table is changed at positions whose weight is less than 

. For every X such that f (X) = 1 in LC, induce 
_ n+1 

f (r X) = 1 on the pure cycles of wt (n-a) >_ 2 , forming 

a heavy cycle (HC) . Note that HC is the reverse complement 

of LC. All of the positions changed so that g (X) = 1 on 

HC are of weights bigger than or equal to 2^- . 

Join LC and HC by using a type RC number, y, which exist 

by Lemma 4, since for n odd the truth table states have 

even (n-1) length. The RC number y has weight . Letting 

f (Y) = 1 then joins LC and HC forming one sequence. By 

Theorem 1 the sequence is RC , and by Lemma 6 the sequence 

is de Bruijn. 

Q.E.D. 
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C. LOWER BOUND ON NUMBER OF RC SEQUENCES 



1 . Theorem 

A weak lower bound on the number of RC de Bruijn 
sequences is established. 

Theorem 7 ; A lower bound on £h| number of de Bruijn 
RC sequences for odd n>3 is A • 2 “ , where A is the 
number of possible different interconnections for pure 
cycles of weight<2^~ as determined by the Best method 
[Ref. 5]. 

Proof: Restricting attention to those cycles of 

wt^S^ (since the other connections are induced) , the 
Best method gives the number A of possible different inter- 
connections fgr^the pure cycles of weight<_2^- >n Jjiince 

2 — 2 

there exist 2 type RC numbers, this gives 2 ways 

of joi|jiijLg A different pairs of cycles for a total of 

2 

A • 2 possibilities. 

Q.E.D. 

2 . Example for n = 7 

The pure cycles of wt<3 and the adjacency matrix 
of A are given in Tables 5.3 and 5.4 respectively. The 
value of the determinant of the adjacency matrix is 38,880 
for any diagonal element. This gives a lower bound of 
38,880 X 8 = 311,040 for the number of de Bruijn RC sequences. 
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0 

1 

3 

5 

7 

9 

11 

13 

19 

21 



TABLE 5.1 



PURE CYCLES FOR N=7 § wt < 3 



2 


4 


8 


16 


32 


64 


6 


12 


24 


48 


96 


65 


10 


20 


40 


80 


33 


66 


14 


28 


56 


112 


97 


67 


18 


36 


72 


17 


34 


68 


22 


44 


88 


49 


98 


69 


26 


52 


104 


81 


35 


70 


38 


76 


25 


50 


100 


73 


42 


84 


41 


82 


37 


74 
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1 

1 

0 

0 

0 

0 

0 

0 

0 

0 



TABLE 5.2 



ADJACENCY MATRIX FOR n =7 



-1 


0 


0 


0 


0 


0 


O ' 


0 


0 


7 


-2 


-2 


-2 


0 


0 


0 


0 


0 


-2 


7 


0 


0 


-2 


-1 


-1 


-1 


0 


-2 


0 


7 


0 


-1 


-1 


-1 


0 


-2 


-2 


0 


0 


7 


0 


-1 


-1 


-2 


-1 


0 


-2 


-1 


0 


5 


0 


0 


0 


0 


0 


-1 


-1 


-1 


0 


3 


0 


0 


0 


0 


-1 


-1 


-1 


0 


0 


3 


0 


0 


0 


-1 


0 


-2 


0 


0 


0 


3 


0 


0 


0 


-2 


-1 


0 


0 


0 


0 


3 
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VI . CONJECTURES 



Although a complete proof for a (c, n) = 0 (mod 4) 
has yet to be given, a technique for grouping RC sequences 
in a 0 (mod 4) fashion will be demonstrated for n = S. 
However, the technique fails to preserve complexity in 
every case. However, for n = 5 the complexitites produced 
were congruent to 0 (mod 2) and when the complement of 
each sequence is added a 0 (mod 4) distribution results. 

It remains, however, that the complexity of the generated 
sequences are not determined apriori. The technique pre- 
sented is interesting in its own right, and it is hoped 
that the interested reader may be able to apply it towards 
a solution of the 0 (mod 4) distribution. 

A. RC SEQUENCE GENERATOR TECHNIQUE 

Table 6.1 shows a listing of four RC sequences and their 
generators. The technique used was to choose in turn each 
of the type RC numbers to have a functional value of 1. 

Note that the complexity of S^ and S^ are both 23, 
while the complexity of S^q and S^ are both 29. At present 
there is no way of knowing which pair of RC sequences will 
have the same complexity, except that it appears C (S^) =C (S^) 
iffall four RC sequences have the same complexity. A listing 
of the generators for each of the 32 pairs of RC sequences 
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TABLE 6.1 







RC GENERATOR 


TECHNIQUE 






C(S.) =23 


C(S 5 )=23 


CCS 10 )=29 


C(S 12 )=29 


n 


GCS.) 


g(s 5 ) 


G(Si 0 ) 


gcs 12 ) 



0 

1 

2 

m 

4 

m 

6 

7 

8 
9 

3D 

11 

~n~ 

13 

14 

15 



1 

1 

0 

1 

0 

1 

1 

1 

1 

0 

0 

0 

1 

1 

1 



1 

1 

0 

0 

1 

_1_ 

1 

1 

1 

1 

0 

0 

0 

1 

1 

1 



1 

1 

0 

0 

1 

0 

1 

1 

1 

1 

_1_ 

0 

0 

1 

1 

1 



1 

1 

0 

0 

1 

0 

1 

1 

1 

1 

0 

0 

J_ 

1 

1 

1 



indicates RC number. 

S 3 = 11111001010001001101110101100000 
S 5 = 11111001011000001101110101000100 
S 10 = 11111001010110000011011101000100 
S 12 = 11111000001101110101100101000100 
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is prepared in Table 6.2. These are listed in order by 
W(G) and then by complexity as far as the groupings allow. 
Keep in mind, that the complement of each RC sequence is not 
included in the table. A close analysis of this table then 
shows a 0 (mod 4) grouping for the number of RC sequences 
of complexity c and span n = 5. 

B. PROBLEMS 

For n = 7, there exist 6 type R numbers (other than 0 

6 64) and 8 type RC numbers in the truth table. This large 

number of possibilities will create considerable problems 

in extending this idea to n = 7 and beyond. 

Using this approach a large scale computer based analysis 

is Essential to any further investigation of RC sequences 

for n>_7 . Considering the number of de Bruijn sequences for 
57 17 

n = 7 to be 2 or 1.42 x 10 an algorithm for generating 

17 

only the RC sequences is crucial, since generating 10 
de Bruijn sequences is not technically feasible. Though 
Theorem 7 establishes the lower bound of 311,040 RC sequences 
for n = 7 , this is a very weak bound and a presumably still 
conservative estimate would be on the order of 10^* actual 
RC sequences for n = 7. 
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TABLE 6.2 



RC de BRUIJN SEQUENCE 0 (Mod 4) 
GROUPING FOR n=5 









G 


(S) 










W[G(S)3 


C(S) 


100 


0 


1 


1 


1001 


0 


0 


1 


101 


7 


25 








0 






1 






7 


25 








1 






IT 






7 


27 




1 




IT 












7 


27 


100 


£) 


0 


0 


1011 


1 


0 


0 


Oil 


7 


27 












IT 




1 




7 


27 




1 












tr 




7 


29 




0 




1 












7 


29 


101 


0 


0 


0 


0010 


0 


1 


1 


Oil 


7 


31 












1 




tr 




7 


31 








1 




0 








7 


31 




1 




IT 












7 


31 


100 


IT 


1 


0 


0010 


0 


0 


1 


111 


7 


31 












1 




IT 




7 


31 








1 




IT 








7 


31 




1 




IT 












7 


31 


101 


0 


0 


0 


1011 


1 


1 


0 


Oil 


9 


31 












0 




1 




9 


51 




1 












IT 




9 


31 




0 




1 












9 


31 


100 


0 


1 


IT 


1011 


1 


0 


0 


111 


9 


31 












0 




1 




9 


31 




1 












IT 




9 


31 




0 




1 












9 


31 


111 


0 


0 


0 


1111 


1 


1 


0 


Oil 


11 


23 












0 




1 




11 


23 




1 












0 




11 


29 




(T 




1 












11 


29 


101 


0 


1 


I 


1011 


0 


1 


0 


111 


11 


25 








0 








1 




11 


25 




1 












IT 




11 


27 




IT 








1 








11 


27 
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APPENDIX A 



ALGORITHM FOR GENERATING LEAST ELEMENT 
FOR EVERY PURE CYCLE 



This algorithm is designed to produce the least element 
of each pure cycle. The least elements are initially in 
decreasing order. A simple rearrangement will group the 
least elements by weight in decreasing order if necessary. 

The algorithm to be presented is an adaptation of the 
0 -operation discussed in a paper by Fredricksen and Maiorana 
[Ref. 10], to generate a lexicographic list of necklaces. 

ALGORITHM: Begin the ft operation with the zero n-tuple,. 

which is the first least element for the zero pure cycle. 

ft Operation: ft (x^ ••• x n ) = (y^ ••• y n ) = Y 

1. Find the largest subscript j such that 



x- =0 and x. = 1 for k>j . 
J K 



2 

3 



Form x^ X 2 



x- ^ 1 where x^ = 1 



Repeat x^ X 2 ••• Xj_-^1 until n numbers are produced, 
i) If n = tj for some integer t then let 



Y=x. 



x j-i lx i' 



x j -1 1 



•N-1 1 



ii) If n>tj then finish with x, •••x . 

J 1 n- t • 



then let Y=x^***x. _^1' 



J 

x j-i lx i 



4 

5 



‘ n 'V 



Y is a least element iff n=tj 
Repeat beginning at step 1 until the final least 



n 



element (1) is reached. 
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TABLE A . 1 



Q ALGORITHM FOR N = 7 

Wt 

0000000 o 

0000001 1 

0000010 * 

0000011 2 

0000100 * 

0000101 2 

0000110 * 

0000111 3 

0001000 * 

0001001 2 

0001010 * 

0001011 3 

0001100 * 

0001101 3 

0001110 * 

0001111 4 

0010010 * 

0010011 3 

0010100 * 

0010101 3 

0010110 * 

0010111 4 

0011001 * 

0011010 * 

0011011 4 

0011100 * 

0011101 4 

0011110 * 

0011111 5 

0101010 * 

0101011 4 

0101101 * 

0101110 * 

0101111 5 

0110110 * 

0110111 5 

0111011 * 

0111101 * 

0111110 * 

0111111 6 

1111111 7 

* Unacceptable 
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Table A.l illustrates the algorithm for n=7. Arranging 
in order by weight, the least element for each pure cycle is 
tabulated in Table A. 2. 



TABLE A. 2 

PURE CYCLE LEAST ELEMENTS (n=7) 



0 


0000000 


0 


1 


0000001 


1 


2 


0000011 


3 




0000101 


5 




0001001 


9 


3 


0000111 


7 




0001011 


11 




0001101 


13 




0010011 


19 




0010101 


21 


4 


0001111 


15 




0010111 


23 




0011011 


27 




0011101 


29 




0101011 


43 


5 


0011111 


31 




0101111 


47 




0110111 


55 


6 


0111111 


63 


7 


1111111 


127 
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